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Abstract 

We consider an inhomogeneous linear Boltzmann equation, with an external confin¬ 
ing potential. The collision operator is a simple relaxation toward a local Maxwellian, 
therefore without diffusion. We prove the exponential time decay toward the global 
Maxwellian, with an explicit rate of decay. The methods are based on hypoelliptic meth¬ 
ods transposed here to get spectral information. They were inspired by former works on 
the Fokker-Planck equation and the main feature of this work is that they are relevant 
although the equation itself has no regularizing properties. 


1 Introduction. 

This article is devoted to the study of the long time behavior of the solutions of the following 
kinetic equation in f of unknown / 

f dtf + v.d^f - d^V(x).dyf = Q(f), , . 

1 /|t=o = /o. ^ ^ 

The right-hand side is a simple linear model for the Boltzmann operator 

Q(f) = liPP-oo - /), p{t, ^) = j ( 2 ) 

where ^oo is the Maxwellian in the velocity direction 

This equation describes a system of large number of particles submitted to an external force 
deriving from a potential V{x), and for which the collision operator in the right-hand side is 
a simple relaxation toward the local Maxwellian pfioo- In particular there is no diffusion. 

We suppose that the derivatives of V of order 2 or more are bounded, and also that 
e~^ G L^, which implies that there is a unique steady state. In this case we say that V is 
a confining potential (anyway the adaptation in the case when e~^ 0 is straightforward, 
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see remark ESI). It can be useful introduce the spatial Maxwellian and global Maxwellian 
defined respectively by 

e-V'W 

^ /e-^Wrfx ’ ^ Poo{x)^ioo{v). 

All steady states in are proportional to the Maxwellian M. In order to study the 

exponential decay we now introduce an additionnal operator, 

Ps} =-^d^{d^ + v)-'^dx{dx + dxV)+ 1. ( 4 ) 

This operator has nice properties in the following weighted space 

g p's, t. G L^idxdv)^ , 

with the natural norm defined by 11/11^2 = / \f\‘^M~^dxdv. Indeed the closure from of 
A2-1 in is maximal accretive (see jl]) and has 0 as single eigenvalue associated with the 
eigenfunction AI. We shall assume the following: 

Operator A^ — 1 has a spectral gap a > 0 in B^. (5) 

Recall that the spectral gap is defined has the infimum of the spectrum except the lowest 
eigenvalue. We mention now some simple cases when it happens. For example when HessR > 
Aid then a = X. It is a special case of the one when |I4'(x)| goes to infinity with x, which 
implies that h? is with compact resolvent in B^ and that © is also satisfied. We refer to 0 
or and reference therein for complementary information about it. 

Now about the collision operator Q, we just mention here that it is mass and positivity 
preserving and ” dissipative” in the sense non-negative in B^ (see PP ). We shall study it more 
carefully later, and refer to remark 14..‘11 here for complements. 

It is easy to verify that —vdx + dxV{x)dy -|- Q is also mass and positivity preserving and 
dissipative, and that its closure in B^ from generates a semi-group of contraction in B^. 
The Cauchy problem o is therefore well posed and it was proven in pQ that under regularity 
assumptions and bounds on the solution / of ©, /(t) tends to M. when t goes to infinity 
faster than any inverse power of t. We now state our main result: 

Theorem 1.1 There exists a constant A > 0 depending only on the second and third order 
derivatives of V, sueh that for all normalized function fo £ B^, we have the following 

||/(T-) - /oo|ls2 < 3 ||/o-/oo|Is2 

here foo = A4, and f is the unique solution of equation m- 

As a direct consequence we also obtain the decrease of the so-called relative entropy: 

Corollary 1.2 Under the hypothesis of the preceding Theorem, and assuming in addition 
that /o > 0, we have 

0 < H{f,foc){t) jj fit) In dxdv < 3||/o||s2 ||/o - /oo|lij 2 


2 




This study is motivated by proving the validity of some new tools and ideas, namely the 
one called hypocoercivity, appeared in a few recent articles in order to prove exponential time- 
decay convergence for some inhomogeneous (mostly linear) kinetic equations such as Fokker- 
Planck inUHIHlIE], Vlasov-Poisson-Fokker-Planck [2j, chains of anharmonic oscillators It 
is well-known that in the homogeneous case the exponential decay can be easily obtained by 
spectral methods if we assume (in the linear case) some coercivity of the collision operator. 
In the inhomogeneous case, the global coercivity is false in general, but can be obtained in 
a modified but norm-equivalent Hilbert space {B^ in this work). This property can serve 
as a definition of hypocoercivity. Several tools can be used for this, essentially inspired by 
hypoelliptic ideas, that’s why this name was introduced very recently in [S]. We mention some 
of them: The use of Kohn’s method to get simultaneously hypocoercivity and hypoellipticity 
0 , 0 , 0 , in the case of the Fokker-Planck operator; the use of analytic dilation and complex 
FBI-Bargmann transform [^, a method of multiplier via pseudodifferential operators or 
functional analysis using harmonic oscillators and Witten Laplacian (linear part of j7j)- 

About the trend to equilibrium, this has been studied for the long time, and we only 
want to quote the entropy dissipating methods introduced by Villani and Desvillettes to prove 
arbitrary and explicit algebraic time decay. It was used for the Fokker-Planck equation j2j, for 
the model studied here PP, and in its main achievement for the full Boltzmann equation 
Now the question naturally arose whether the exponential decay, obtained via hypoelliptic 
tools is also true for non-hypoelliptic operators. In this work we choose a simple example 
of collision operator which has no regularity property. It appears that Lie techniques (those 
also in the core of the hypoelliptic theory) also give sufficient information on the spectrum 
and in particular the spectral gap (in some modified space) to get hypocoercivity an then 
exponential decay. We therefore hope this techniques to be applied in the future to other 
inhomogeneous kinetic equations with linear or non-linear collision operators with or without 
regularity properties (see the review [0] for examples). 

Eventually our result answer a question raised by Caceres, Carillo and Goudon in PP 
about the applicability of hypoelliptic techniques of [H] to obtain explicit exponential decay 
of the model studied here. 


Contents 


1 

Introduction. 

2 

Notation.s and functional analvsisl 


3 

Hypocoercivity for operator K. 


4 

Proof of the Theorem and comments! 


Bibliography 


1 

3 

5 

9 

10 


2 Notations and functional analysis 

For a Hilbert-space type analysis, it is more convenient to work directly on the equation 
obtained after conjugating with Therefore we pose u = which is now supposed 
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to vary in and the equation satisfied by u is then 


r dtu + Ku = 0 

\ u\t=o = uo, ^ ^ 

where we introduced the notations 

' K = Xo + 7(Id-ni), 

, Xo = v.dx - dxV(x).dy, ('y) 

niu(x,u) = (^f u(x,v')/ii!^^(v')dv''j for u G L^. 

It is immediate to see that Hi is an orthogonal projector onto 

-El = |u s.t. exists p G L?{dx) with u = 


which is a closed subspace of L^. We therefore directly get that in the new formulation in 
-L^, the collision operator, which is now —(Id — Hi) is dissipative. Of course we also have 


so that the square root of Maxwellian is in the kernel of the (new) collision operator. It is 
clear that K and its adjoint K* = —Xq + 7 (Id — Hi) are well defined in in 

and as (non closed) operators in L^(M^'^) with domain D{K) = D{K*) = 5(M^'^). We denote 
by the same later their maximal closure in L^. For j = 1,... ,d, we introduce the differential 
operators aj, the annihilation operator bj‘. 

= 7 ^'^^ {dxj + dx^V{x)/2) bj = {dy. + Vj/2) , 
and their formal adjoints 

a* = 7 ^/^ {-dxj + dx^V{x)/2) b* = 7 ^/^ {-dy. + Vj/2) . 

For the sake of conciseness, the letters a and b denote the vectors 



/ Oi \ 

( \ 

0 = 

: 

\ad / 

UJ 


while a* and b* are the forms 


a* = {al,...,a*,) b* = {b*,,... ,b*,). 

Up to the factor 7 , the non negative operator — 1 is nothing but the sum of the Witten 
Laplacian on 0-forms 7 “^a*a = —Xx + \dxV{x)^ /4 —AU(x)/2 and of the harmonic oscillator 
'y~^b*b = — -|- u^/4 — d/2. Under the hypothesis on U, — 1 is a nonnegative, is 

a core and A^ is well-defined for r G M. Moreover the kernel of A^ — 1 is Eo = CM 1/2 _ Yot 
the following we shall denote IIo the orthogonal projector onto Eq, 


Uqu = 

which also extends to n G 5'(M'^). Let us now recall some relations involving the operators 
a, b and Xq, that can be found for example in section 1 of We mention here that since 
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the operators are continuous in S and S' there is no problem of defining their commutators. 
First recall the Canonical Commutation Relations (CCR) for b and their counterparts for a 

[bj ,bk] = [bj , 6fc] = 0 [bj, bl] = 'ySjk [aj, o^] = [a^, aj] = 0 [aj, al] = V. (8) 

It is also clear that the a’s and a*’s commute with the b’s and b*’s. The main remark is that 
the Uj’s, Qj’s are in the Lie algebra generated by the 5j’s, b*’s and the vector field Xq: 


lbj,Xoj = aj [h*,Xo] = a*. ( 9 ) 

Similarly, the 6 j’s and 6 *’s can be derived from the a^’s, a*’s and Xq 

d d 

[aj,XQ] = -Yl ^o] = • (10) 

k=l k=l 

The relations and m are summarized by 

[b,Xo] = a , [b*,Xo] = a* , [a,Xo] = -HessRb and [a*,Xo] = - 6 *HessR. (11) 

By combination we have the useful formulas: 

[A^, Xo] = - 6 *(HessR - Id)a - a*(HessR - Id) 6 , (12) 

b*{a*a) = {a*a)b* a*{a*a) = {a*a)a* - ya^HessR^ (13) 

and a*{b*b) = ib*b)a* b*{b*b) = ( 6 * 6 ) 6 * - 76 *, (14) 


and their adjoint relations hold as equality of continuous operators in 5(M^'^) and S'{M.‘^'^). 


3 Hypocoercivity for operator K. 

In this section we continue to work with operator K defined in the preceding section. We 
shall prove that operator K is hypocoercitive, i.e. coercitive in with a modified scalar 
product. For this we introduce an additional operator 

L = A-^a*b = A-2(^ a*bj) 
j 

We shall see later that this operator is explicitly bounded in terms of the second and third 
derivatives of V. 

Proposition 3.1 Let a be defined in 0). Then there exists e, A > 0 such that for all 
S 9 uTM1/2 we have 

o? o 

Re (iLtt, {Id + e{L + L*))u) > — ||ri|| , 

./I 

where A can he chosen to depend explicitly on 7 and the second and third derivatives of V, 
and \\eL\\ < 1. 
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Proof. Let us take u £ L'^ and e > 0. We write 


Re {Ku, (Id + e(L + L*))u) 

= Re 7 ((Id - ni)u, (Id + e{L + L*))u) + Re (Xqu, (Id + e{L + L*))u) 

= 7 II(Id - ni)uf + e 7 Re ((Id - ni)u, (L + L*)u) + eRe {Xqu, +{L + L*)u) 

= / + // + III, 

where in the last term we used the fact that Xq is skew-adjoint. We first study the two first 
terms. Using the Cauchy-Schwartz inequality we can write 

I + II>^\\{Id- Rijnf - ||(i + i*)^||2 > 1 ||(J^ _ n,)^||2 _ ^2^ ||^||2 ||^||2 ^ 

Now we study more carefully the third one 

III = eRe {Xqu, (L + L*)u) = eRe ([L, Xo]u, u) , (16) 

again since Xq is skewadjoint. Recalling that L = A~‘^a*b we can write using the equalities 

m 

[L,Xo] = [A-\*b,Xo] 

= [A-^,Xo]a*b + A-WXo]b + A-^a*[b,Xo] (17) 

= -A-^[A‘^,Xo]A-^a% - A-H*RessVb + A-^a*a. 

Here we used the fact that for A, B, and B~^ continuous on S and S' we have [A,B~^] = 
— B~^[A, B]B~^. Let us denote 

A = -A-^[A^,Xo]A-‘^a*b - A-‘^b*RessVb 

We postpone to the end of this section the proof of following lemma concerning operators A 
and L. 

Lemma 3.2 Operator A and L are bounded on L^. Besides their norms of can be explicitly 
bounded in terms of j and the second and third order derivatives ofV. 

Now it is clear that denoting for u G L^}, operator b considered as an operator 

from if into II^~^ satisfies Ei C Ker(6). In fact this comes from the fact that b is then 
the annihilation operator in the velocity variable. Since it appears only in the right in the 
expression of A we therefore get 

El C Ker(.4). 

As a consequence we can write A = A{I — Hi) and therefore 

[L, Xo] = .A(Id - Hi) + A-^a*a 

Putting this in the expression of the term III appearing in m yields 

III = eRe (^(Id — ni)u, u) -)- eRe [A~^a*au, u) 

2 MS') 

> — Y ||(Id — ni)u||^-ll-^ll^ ll^ll^ + Ale (A~‘^a*au, u) . 

4 7 ' 
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Now it is also clear that A^, a*a and Pi commute together, and we can therefore write for 
the second term appearing in the last inequality 

eRe (A“^a*att, u) = eRe (yh~^a*alliu,u) + eRe (A“^a*a(l — ni)u, u) 

= eRe (A“^a*aniu, IIiu) + eRe (A“^a*a(l — ni)n, (1 — ni)u) (19) 
> eRe (A“^a*anin, Ifiu) — e || (1 — ni)n|p . 

for the last inequality we simply used the fact that a* a < A^, which implies easily that the 
norm of oA“^ is bounded by 1 . 

Now we can take into account the spectral gap property of A^ together with the fact that 
u is supposed to orthogonal to We write it as a lemma 

Lemma 3.3 Recalling that u A we have Re (A“^a*anin, Ilin) > ll^ll^ 

Proof. We first notice that A^ and a*a leave Pi invariant, and that 

= 1 + a*a 

is essentially the Witten Laplacian in the spacial direction. We define now r to be the spectral 
gap for a*a. Now for w G 5(M^) such that w A p\!^ we have 

{a*aw, re)i 2 (]Rd) > r ||'w||i 2 (iRd) . 

Since {a*a A A also we get 

{a*a{a*a + {a*a + ■ 

Now the following inequalities are clear; 

T a a 

- > -> - . 

lAr lAo lAy 

Indeed t > a from the definition of a and a < 7 because of the harmonic part of A^ — 1. Now 
since Hitt G Pi n Pq- we get the result of the lemma from the preceding study by applying it 
to the function defined for a.e. u by x 1 —> Iiiu{x,v). □ 

End of the proof of Proposition 1.3,1 L Now we can put the result of the lemma in (|19ll 
and we get 

eRe [A~^a*au,u) A — 11 ^ 11 * 11 ^ — e || (1 — ni)u||^ . (20) 

We obtain the following lower bound for the term III from m 

III > ||(Id - ni)uf - - Pf ||uf A llRiuf - e 11(1 - ni)nf . (21) 

4 7 1A 7 

Eventually putting together the estimate on lAlI and III we get 

Re {Ku, (Id Ae(L AL*))n) 

>|||(/<i-ni)ti||"-e 7 ||L||"||uf 

-1 ||(id - ni)uf - d Mil" Hull" + iln.uli" - E 11(1 - nouli". 

4 7 0 ^ A a 

> I ||(Id - ni)uf A pi^xf - ( 7 "^ Pf A 7 Pf) pf 
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by taking e < 7 / 8 . Now we use the fact that Hi is an orthogonal projector and that s < 7 / 8 : 

Re {Ku, (Id + e{L + L*))u) > Pf + 7 Pf)) hf ■ 

For e/a sufficiently small, but depending only on 7 and the second and third order derivatives 
of V via lemma we obtain 

a'^ 9 

Re (Ku, (Id + e{L + L*))u) > — ||tt|| 

where A satisfies the hypothesis of the Proposition. The proof is then complete. □ 


Proof of lemma D Recall that L = A ‘^a*b. Now for A we use the expression of 
[A^^o] in m and we get 

A =A-‘^b*(Ress{V) - ld)aA-^a*b + A-^a*(Ress{V) - ld)bA-^a*b - A-^b*Ress{V)b. ( 22 ) 


We therefore see that it is sufficient to prove that for any dx d real matrix M{x) depending 
only on x, bounded and with first derivative bounded, the following operators 

A-H*M{x)a, A-H*M{x)b, A-^a*M{x)b 

are bounded as operators on We give the proof for the first one since for the remaining 
ones, the proof is similar and easier. We shall prove the result for its adjoint a*M{x)bA~‘^. 
We write for u G 5, 


*M{x)bA ^u|| < 


j,k 


<'^\\Mj,kix)a*bkA ^mII + ^7^/2 ||( 9 a,^.Af,-fc)(a;) 6 fcA ^23) 


< 


+ 7 ^/^ \\d^M\\^^j ^ {\\a*jbkA ^ri|| + || 6 fcA 


- 2 „ 


where we used the fact that [a*,Mj^k] = —dxjMj^k- Now this is straightforward to check that 
|| 6 fcA“^n|| < ||u|| since b*b < A^ and 1 < A^. For the other term 
inequality of (P|) , we write 


a*bkA “^u 


in the last 


||a*6fcA = [aja*bkA '^u,bk^ 

< {a*ajbkA~'^u,bkA~‘^u) + (^j{dl.V)bkA~‘^u,bkA~'^u^ 

< (A^6fcA“^n, 6fcA“^n) + 7 ||HessP||^oo [bkA~‘^u,bkA~‘^u) 

< (A^6fcA“^n, 6fcA“^n) + 7 ||HessI/||^oo {bkA~‘^u,bkA~‘^u) 


Now using the fact that [A^, bk] = —jbk we can continue the preceding series of inequalities: 

< [A^bkA~‘^u,bkA~‘^u) + 7 ||HessP||^oo {bkA~‘^u,bkA~‘^u) 

< {bkU,bkA~‘^u) +7(||HessI/||^oo + 1) {bkA~‘^u,bkA~‘^u) 

< 7(||HessI/||^oo + 2) ||ti|p 
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in the last inequality 

the second and third 

□ 


a*bkA 


-2 


U 


again since b’^b^ < and 1 < A^. Therefore the term 

of (f^ is also bounded by C'HttH, where C depends only on 7 anc 
derivatives of V. The proof of lemma is therefore complete. 


4 Proof of the Theorem and comments 

We go on studying operator K defined in the preceding sections. We first quote an easy 
result from [Jj relying hypocoercivity to exponential decay 

Lemma 4.1 (lemma A6 in [ 7 |) Let JC be the infinitesimal generator of a semigroup of 
contraction on a Hilbert space H and suppose that there exist a constant (5 > 0 and a bounded 
operator C with norm bounded by Cc > 1 such that 

VuET»(/C), (5||uf < Re(/C7r,n)+Re(/Cn,(£ + £*)u) (25) 

6t 

then for all uo £ H and t > 0 we have ||e“*^rto|| < 3e ||riQ|| 

Proof of Theorem 11.11 we first consider operator K defined in the preceding sec¬ 
tions.The result for K in the Hilbert space is then a direct consequence of Propo¬ 

sition o Indeed it suffices to apply the preceding lemma when considering K in place of 
JC, eL in place of C, Cc = 1 and /A in place of 6. The result of the theorem is then 
a simple transcription in terms of / = and / — /oo = — Hqu) (of course 

{u — Holt) T which is an other way to say that fj(f — foo)dxdv = 0). The proof is 

complete. □ 


Proof of Corollary II.21 This is then a direct consequence of the main Theorem, and the 
proof follows exactly the one in j^l corollary 0.2]: Let /o be a L^-normalized function which 
belongs to and let /(t) be the solution of (pQ), (it stays non-negative for all time). 

The non-negativity of the relative entropy is clear from the fact that / stays normalized. 
For the other side, we write for t > 0, with /oo = Af: 

H{f{t)\M)= 11 /(t)ln(^^^ dxdv = jj dxdv. 

Applying first ln(x) < x — 1 and then the Cauchy-Schwarz inequality for t > 0, we get 


H{f{t)\M) < 


fit) 

AIV2 


< 


fit) 


A1V2 


fit) 


m 

M 


— 1 I dxdv 


AII/2 




(26) 


= \\fit)\\B^ \\fit)-M\\i^2 . 


Hence the exponential decay of the relative entropy is a consequence of Theorem ll.il . 


□ 


We end this work by making some remarks about the main result: 

Remark 4.2 In the case e~^ 0 L^, Theorem II. II remains true when replacing /oo by 0. 
A careful study of the proof shows that it is exactly the same with the following adapta¬ 
tions: Pooix) = is not anymore in L^{dx) and neither does Ai[x,v) = Pooix)pooiv) 
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in L^{dxdv). Then a is the bottom of the spectrum of operator — 1, and in the proof 
Eo = {0} and the corresponding projector is IIo = 0. In this case Theorem 11.11 has to be 
understood has a vanishing. 

Remark 4.3 A natural question is to understand the common features between for example 
the Fokker-Planck operator (and also the chains of anaharmonic oscillators) and the one 
studied here. Let us see this after conjugating by the square root of the Maxwellian and only 
in for the velocity variable (homogeneous case). The collision operators are respectively 

Qfp = b*b (Fokker-Planck) , Quf, = Id — Hi (Linear Inhom. Boltzmann). 

Here Hi is the orthogonal projection on the space spanned by the Maxwellian in the velocity 
variable /ioo • A simple remark can be made using the Hermite decomposition of functions 
(in the velocity variable), which we denote We know that the annihilation operators h 
and its conjugate b* (creation operator) have both a nice description, since for all /c > 1, 

bHk = VkHk-u b*Hk-i = y/kHk 

and bHo = 0 (recall that Hq = /Too )• Therefore b and b* can be represented by respectively 
an upperdiagonal and a subdiagonal infinite matrices with coefficients Immediately we 
get the expression of the harmonic oscillator b*b, for which the Hermite polynomials are an 
Hilbertian base. Now using this decomposition we can build operators c and its adjoint c* 
by imposing 

cHk = Hk-i, = Hk 

and cHq = 0. Then c and c* have the same representation as matrices than b and b* 
respectively (note that they are bounded). It is immediate to check that 

c*c = Id - Hi = Qiih. 

As a conclusion, and transferring this in the inhomogeneous case, we see that the Fokker- 
Planck and the linear inhomogeneous Boltzmann models have the same structure, explaining 
(a bit) the efficiency of hypoelliptic methods in the last case. 
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